Small pitch-angle magnetobremsstrahlung in inhomogeneous curved magnetic fields 



S.R. Kelner[| and F.A. Aharonian| 

The character of radiation of relativistic charged particles in strong magnetic fields largely depends 
on the disposition of particle trajectories relative to the field lines. The motion of particles with 
trajectories close to the curved magnetic lines is usually referred to the so-called curvature radiation. 
The latter is treated within the formalism of synchrotron radiation by replacing the particle Larmor 
radius with the curvature radius of the field lines. However, even at small pitch angles, the curvatures 
of the particle trajectory and the field line may differ significantly. Moreover, as we show in this 
paper the trajectory curvature varies with time, i.e. the process has a stochastic character. Therefore 
for calculations of observable characteristics of radiation by an ensemble of particles, the radiation 
intensities should be averaged over time. In this paper, for determination of particle trajectories 
we use the Hamiltonian formalism, and show that close to curved magnetic lines, for the given 
configuration of the magnetic field, the initial point and particle energy, always exist a smooth 
trajectory without fast oscillations of the curvature radius. This is the trajectory which is responsible 
for the curvature radiation. The realization of this regime requires the initial particle velocity to 
be directed strictly along the smooth trajectory. If initially the velocity vector is adjacent to the 
field line or close to it, a transition of the particle to the smooth trajectory on a later stage cannot 
be excluded, however under certain conditions this can happen only after the particle has already 
radiated out all its energy. Even for an initial narrow angular distribution of the particle beam, 
the energy spectrum of the resulting radiation appears significantly harder compared to the energy 
spectrum of the curvature radiation. This result might have direct relation to the recent spectral 
measurements of gamma-radiation of pulsars by the Fermi Gamma-ray Space Telescope. 

PACS numbers: 41.60.Ap, 41.60.-m, 95.30.Gv 



I. INTRODUCTION 



Eq.(m) follows that 



The motion of a relativistic charged particle in a strong 
magnetic field is accompanied by a radiation which 
in general terms can be called magnetobremsstrahlung. 
While moving in a homogeneous field, the particle energy 
is gradually radiated away, however, its velocity compo- 
nent I'll parallel to the magnetic field remains unchanged. 
Indeed, this statement is obvious for = 0, therefore in 
an arbitrary coordinate system which moves along (or 
opposite) the magnetic field, we have = const. There- 
fore, 
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where and are the parallel and perpendicular com- 
ponents of the momentum, respectively. This equation 
allows us to find the remaining energy of the particle, 
after the radiative damping of the perpendicular compo- 
nent of motion. Let us denote the initial and final values 
of the parallel and perpendicular components of the mo- 
mentum as p^, and p'^ = 0, respectively. From 
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If the initial perpendicular momentum is nonrelativistic, 
Pj_ <C mc, then pj| ^ , i.e. the damping of the perpen- 
dicular motion does not have an impact on the particle 
energy However, for mc <C <^ 7 we have pj| <C Py , 
i.e. during the damping of the perpendicular motion the 
particle loses almost all its energy At ^ mc, the mo- 
menta py and p^ vary similarly with time since the ratio 
Pj_/P\\ ~ const (see Eq.(^. The constant ratio of p^/py 
has the following explanation. The ultrarelativistic parti- 
cle emits photons within the angle < I/7 in the direction 
of its momentum p, therefore in the case of 1 /7 <C /py , 
one may assume that the recoil momentum is directed 
opposite to the vector p. 

The radiation of charged particles in strong magnetic 
fields is determined by the character of their motion, in 
particular by the disposition of particle trajectories rela- 
tive to the field lines. As we show in this paper, in the 
close vicinity of curved field lines should exist a special 
trajectory the particle motion along which is not accom- 
panied by perpendicular oscillation. We call this trajec- 
tory as smooth trajectory. In the case of homogeneous 
magnetic field the smooth trajectory is a straight line 
parallel to the field line. The initial short interval of 
the smooth trajectory in the curved magnetic field can 
be treated as a circle with a radius close to the curvature 
radius of the field line. An ultrarelativistic particle moves 
steadily along this circle, therefore for calculations of the 
accompanying radiation one can apply the well known 
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formulae, and the radiation itself can be considered as a 
curvature radiation. 

In the case of large pitch-angles in the curved inhomo- 
geneous magnetic field, the inhomogeneity does not affect 
the radiation if the Larmor radius is small compared to 
the curvature of the field lines, and to the characteristic 
distance scale on which the field is noticeably changed. 
This case corresponds to the regime of synchrotron ra- 
diation. All formulae that describe the synchrotron and 
curvature radiation regimes coincide if we present them 
in the terms of the trajectory curvature radius Rc. 

For small pitch angles the picture can be essentially 
different. Let us assume that initially the velocity vec- 
tor is parallel to the field line. Then, the accelera- 
tion = 0, correspondingly the curvature radius[| 
^c|^^Q = CO, thus no radiation take place at t = 0. In 
the curved magnetic field, a rectilinear uniform motion 
is not possible, therefore an acceleration and radiation of 
the particle are unavoidable, and, what is essential for 
the correct treatment of the problem, the acceleration 
and the curvature radius appear to be time-dependent. 
In this paper we show that this is true for all trajectories 
located close to the smooth trajectory. This implies that 
the radiation also varies with time, thus for calculations 
of measurable characteristics of radiation one should av- 
erage the relevant distributions over time. 

It is generally believed that for the initially small pitch 
angles the charged particles rapidly radiate away the frac- 
tion of their energy related to the component of the mo- 
mentum perpendicular to the field line, i.e. the perpen- 
dicular motion quickly disappears, and particles move 
along the field lines (see e.g. Ref. [Q). However, here 
one should distinguish two different cases. Let and 
be the particle momentum components parallel and per- 
pendicular to the smooth trajectory, respectively. Let 
us assume that mc and ^ p^. If <C mc, 

then the perpendicular component of motion is quickly 
damped; the particle will continue to move along the 
smooth trajectory and thus emit curvature radiation. 
However, if mc^ then as in the case of the uni- 

form magnetic field, and p^ will be reduced with the 
same time-dependency. The calculations of character- 
istics of the resulting radiation, which is different from 
both synchrotron and curvature radiation components, 
is the prime objective of this paper. We call this regime 
of radiation as small pitch angle magnetobremsstrahlung. 



magnetic field has the same symmetry as the field of 
an infinitely long straight wire. Then, in the cylindri- 
cal system of coordinates, (r, ^, z), the field has only an 
azimuthal component; it depends only on the distance to 
the axis of symmetry, r, but not on z and (j). We denote 
by (cr, 60, Cz) a moving system of unit vectors linked to 
the point (r, ^, z). The magnetic field can be expressed 
in the form B = B{r) e^. The field lines are circles of ra- 
dius r. It appears that particles in such a field can move 
along trajectories close to the circular field lines. 

The equation that describes the particle motion in the 
magnetic field is 



{v X B). 



(3) 



Let us find a solution for which r = ro = const, and the 
velocity 



V(f)e(j: 



(4) 



Note that the components and Vz are constants. The 
vector 60 is uniformly rotating with an angular frequency 
^ = V(f)/rQ^ therefore the derivative de(f)/dt = —v^Cr/r, 
and the vector = const. Substituting Eqs. (^) into 
(0), we find 



^0 e 

mj— Br = - Vz-D Br • 

ro c 

From here follows that 

Vz mcv^j 
Vcfy eBro 



(5) 



(6) 



Writing the velocity components in the form Vz = v sin a, 
V(f) = vcosa (the angle starts from the "equator", i.e. 
— 7r/2 < a < 7r/2), and expressing the velocity v via 
Lorentz factor, Eq.(^) can be presented in the following 
form 



sma 



mc 7 
cos^ a eB ro 



(7) 



Apparently, for any given 7 this equation defines the an- 
gle a. Eq.(||) indeed has solutions in the form of Eq.(^) 
which describe the helicoidal motion depending of pa- 
rameters 7 and ro- Note that motion of particles with 
strictly circular trajectories is not possible. 

The condition for a small step in the helix, \a\ <C 1, 
requires 



II. THE SMOOTH AND CLOSE TO THEM 
TRAJECTORIES 

We start with the analysis of motion of ultrarelativis- 
tic particle in a simple model which allows exact solu- 
tions. Namely, below we consider the case when the 



^ Note that at the motion with a constant absolute speed value, 
the curvature radius is equal, by definition, to Rc = v'^/\a\. 



K 7 < 



\e\Bro 



Note that typical for pulsars values of B and ro readily 
satisfy this condition. Then, 



eBro ' 



lal <Cl, 



(9) 



thus the particle trajectory can be interpreted as a mo- 
tion along the field line with a simultaneous drift in 
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the perpendicular direction. The positive and negative 
charges drift in opposite directions. Under the condition 
of Eq.(||), the curvature radius Rc of a particle located 
at the distance ro, does not depend on the particle en- 
ergy and practically coincides with the curvature radius 
of the field line ro- Note that this conclusion is correct 
for any configuration of the magnetic field. Indeed, the 
short segments of the field line can be treated as a circle, 
therefore when satisfying the condition of Eq. (||) , always 
exist trajectories which (or at least their initial segments) 
are very close (but not identical) to the field line. In the 
Introduction, such trajectories have been called smooth 
trajectories. Also as indicated above, the magnetobrem- 
strahlung of a particle moving along a smooth trajectory 
should be considered as the strict definition of the cur- 
vature radiation. Its intensity is defined as 



2e^c 7^ 



2e^c 7^ 



(10) 



When the particle moves close to but not exactly along 
the smooth trajectories, the radiation intensity can de- 
viate from Eq.(lO). Therefore one should investigate the 
trajectories in the vicinity of the smooth trajectories. 

For this purpose it is convenient to use the Hamilto- 
nian formalism. The vector potential can be taken in 



the form Ax = 0, Ay 



0, A, 



A{r) 



therefore the 

azimuthal component of the field B{r) = —dA/dr. In 
the cylindrical coordinates the Hamilton function can be 
written as 



H = c\ P^^^^iP,-- A(r))2 + m2c2 , (11) 



where P^, and P^ are the generalized momenta corre- 
sponding to the coordinates r, (j) and z. Because of the 
azimuthal symmetry of the magnetic field and its homo- 
geneity along the z axis, (j) and z are cyclic coordinates. 
Therefore P^ and P^, as well as the energy, are integrals 
of motion. The presence of three integrals of motion al- 
lows us to reduce the problem to the one-dimensional 
case, and thus comprehensively study the qualitative be- 
havior of the solutions. Below we will limit the treatment 
of the problem by the case of ultrarelativistic particles 
and small pitch- angles. 

For given P^ and P^, the radial motion can be consid- 
ered, as it follows from Eq.(pT[), as a motion in the field 
with the following effective potential: 



C/eff(r) = 



(12) 



The first and second derivatives of this potential are 

U^Ar) = ^ + -B{rm - - A{r)) , (13) 

r'^ c c 



and 



2e^ 



U:Ar) = -f + -B'(r){P,— A(r)) + -^B\r). (14) 



Assuming that at r = ro the derivative U^f^{ro) = 0, then 
at this point 



f/:ff(ro) = 



2Pl 



(15) 



B{ro) 



B\ro) . (16) 



If B decreases with r as a power-law, B ~ l/r"^, then 
for J < 3 the both terms in Eq.(^6|) are positive, therefore 
/7eff has a minimum at ro- Furthermore, assuming that 
the condition of Eq.(^) is fulfilled, in Eq.(^6|) we can keep 
only the last term. 

For integration of the equations of motions we apply 
the standard method which is used in the theory of small 
oscillations. In the vicinity of the minimum point the 
function t/eff(^) is approximated by a parabola: 



U,E{r)=P^^m'ujt{r-r^) 



where 



2m? 



eB{ro 



mc 



Pi = C/eff(ro) 



(17) 



(18) 



(19) 



In these denotations, the Hamilton function can be pre- 
sented in the form 



H = c^JP^ + m2a;2(r - ro)^ + P§ + m^c^ , (20) 

and correspondingly the equations of motion are 

dH c^Pr • dH m^c^ujl , 
dPr E dr E 

(21) 

where E is the particle energy. The solution of these 
equations gives 

r = To — pcos , Pr = mpujc sm , (22) 

7 7 

where 7 = E/{mc^) is the Lorentz factor of the particle, 
p is an arbitrary constant^ related to the particle energy 
through the relation 



E = c\ Pq + m?p'^uj1 + m?c^ . 



(23) 



^ Note that here the origin of the coordinate system is arbitrary, 
therefore instead of t one can use t — to, where to is the second 
arbitrary constant. Below we omit the arbitrary constants which 
can be added to 6 and z. 
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The condition of applicability of Eq. ( |22| ) is the smallness 
of p compared to ro as well as to the characteristic linear 
scale on which the magnetic field is changed significantly. 

Now we should find the time-dependencies of other co- 
ordinates. In this regard we note that 



dH 

dP6 



Ps 1 



m7 

Using Eq.dl), we find 



m7 Tq 



2(r-ro) 



Pch 1 



m7 Tq 

Now substituting in 
dH 



^ 1 + — 



2p ujct 
cos 

ro 7 

2p7 . uJct 
sm 

^c^o 7 



(24) 

(25) 
(26) 

(27) 



the following approximate presentation for A{r), 
A{r) = A{rQ + (r - ro)) « A(ro) - B(ro) ■ (r - ro) , (28) 



and using Eqs. (|15[ ) and (P2D, we find 
Pi 



COS 

7 7 



t — p sm , 

7 



(29) 



(30) 



In order to make the analytical expressions more con- 
venient to work with, it is useful to use instead of and 
p the following parameters /^y and I3j_ which are deter- 
mined as 



/?ii = 



mcjro 



P. = ^, (31) 

7C 



and denote 



ro Wc 



(32) 



Since is the projection of the particle momentum 
on the axis then /^y is the component of the velocity 
parallel to the magnetic field (in units of c), and is the 
perpendicular components of the velocity. Finally is 
the drift velocity. 

For the new variables we find 

r/c = /3^ sin r, i/c = — cost , (33) 



(34) 



where r = uJctM] azimuthal component of velocity is 

«^ = r^ = c/3|| M cosrj . (35) 



For = 0, this solution coincides with the obtained 
above exact solution that describes the motion along a 
smooth trajectory (in this case, along the helicoidal line). 
The same trajectory can be found if we average Eqs.(^)- 
( ^5| ) over time. Therefore this solution can be interpreted 
as a motion along the spiral around the smooth trajec- 
tory. At 7^ 0, the solution is approximate; it is correct 
when the following condition is fulfilled: 



(36) 



The relation between and can be arbitrary. Since 
^1, the second of the above conditions is equivalent 
to Eq.(||). The same requirement can be formulated in 
terms smallness of the Larmor radius compared to the 
curvature radius of the field line. 

Note that the inclusion of the process of the gradi- 
ent drift into consideration hardly can be justified since 
it would exceed the accuracy of the approach. Indeed, 
the gradient drift would appear in the formulae if one 
adds in Eq. (p8|) the next (quadratic) term of expansion, 
while our study is limited for the linear terms. Therefore 
throughout the paper we will take into account only the 
curvature drift. 

For the chosen new variables, the velocity and its 
square are 



V — rCr + r(j)ed 



ZCz 



and 



(37) 



(38) 



Note that v'^ does not depend on t. In order to find the 
acceleration, we use the following equation 

a = (v X B) = — (re^ - zCr) , (39) 
mc7 7 

where it is taken into account that the field has only an 
azimuthal component. The square of the acceleration 

= + ,40) 

where = c^/3y/ro. This equation can be derived also 
through the consideration of kinematics. Indeed, 

a = r = (f — r(j)'^)er + (2 + r(j))eci) + zBz . (41) 

From here, due to the smallness of the ratios Pj_/P^^ and 
Pd/P\\^ we obtain an equation identical to Eq.(|40|). 

At = we have = = const. This is the case 
when the particle moves along the smooth trajectory. For 
P± = Pd^ the acceleration = 2ao (1 — cost). For the 
moments of time, when the acceleration becomes zero, 
the particle velocity is parallel to the field line. In the 
limit of /3j_ the acceleration does not depend on 

the radius of the field curvature and on time: 



a = an 



0i 



(42) 
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It is well known (see e.g. Ref. that this equation 
describes also the acceleration of the particle in the ho- 
mogeneous magnetic field. For applicability of Eq.(|2|) it 
is sufficient to fulfill only the condition ^ while 
the relation between and /^y can be arbitrary. 

For an ultrarelativistic particle, the curvature radius of 
trajectory and the acceleration are linked with a simple 
relation \a\ = (? jRc^ which at v = c is, in fact, the 
definition of the curvature radius. Thus, the radius of 
the curvature of trajectory is 

i?eW = ro(^l-^COST+||) ' . (43) 

Here /^y is set to be equal to 1 given that 1 — /^y <C 1. 

Let us consider now the motion of the charged particle 
in the field of a magnetic dipole (the so-called Stormer 
problem). Such a system appears to be (unlike the previ- 
ous case) non-integrable (see e.g. Ref. [||), and generally 
it is impossible to conduct even a qualitative analysis of 
the motion. However, when the condition of Eq.(^) is 
fulfilled, the system has smooth trajectories. The tra- 
jectories in the close vicinity of the latter present lines 
spiraling around the smooth trajectories, i.e. basically 
we deal with the same picture as discussed above for the 
straight wire. This is confirmed also by our numerical 
calculations. 

In Fig. 1^ we show the initial part of the trajectory in 
the field of magnetic dipole directed towards the z-axis. 
At the initial moment of time the particle is located in 
the point r(0) = (sin'i^, 0, cos'i^) where the coordinates 
are measured in units of the radius of the object (star), 
i?*, the angle = 20°, the initial velocity is directed 
along the field line, and the parameter ^ = 10~^ (see 
Appendix ^). The amplitude of oscillations appears of 
the same order of magnitude, i.e. 10~^, therefore in the 
left panel of Fig. |l| the oscillations are not visible. Fig. || 
shows the dependence of acceleration on time. The pe- 
riod of oscillations is equal to Tq = 27rmcj/{eBo)^ where 
^0 is the magnetic field on the surface of the star. The 
decrease of the amplitude and the increase of the period 
is due to the reduction of the field along the trajectory. 
The change of the Lorentz factor 7 due to the radiation 
losses is not taken into account. The trajectory patters 
for larger time intervals are discussed in Appendix. 

It should be noted that the necessity of calculations of 
characteristics of (synchro-bremsstrahlung) radiation in 
the context of proper treatment of trajectories of charged 
particles in curved magnetic fields, has been correctly in- 
dicated also in Ref. [Q. However, our analysis shows 
that some of the basic results and conclusions of that pa- 
per appear to be not correct. Therefore for clarification 
of some issues of conceptual character, below we com- 
pare the results of the present study with the relevant 
result from Ref. Q. For convenience, the equations from 
Ref. m will be referred by their original numbers with 
an additional indication of Q]. 

Let us compare, in particular, the results of calcula- 
tions regarding the square of the acceleration, a^, in the 



z y * 2.7*10^ 




0.35 0.36 0.35 0.36 



X X 

FIG. 1. The initial segment of the trajectory of charged par- 
ticle in the field of the magnetic dipole. The initial point is 
located on the surface of the star at the angle 20° relative 
to the ^-axis, the initial velocity is directed along the field 
line. The y-axis in the right panes is rescaled (multiplied by 
a factor of 2.7 x 10^ in order to make visible the oscillations. 



[arb. units] 




FIG. 2. The time-dependence of acceleration during the mo- 
tion of a charged particle along the trajectory presented in 

Fig.g 



curved magnetic field. For the denotations used in this 
paper, the equation (5,Q) can be written in the form 

a^=«S(l + |i)\ (44) 

This equation is not correct; in particular it lacks time- 
dependence. Apparently, the error in Ref. [|| has ap- 
peared because of the definition of acceleration for which 
the authors used a wrong formula given by equation 
(1,1^). The reference to the monograph of Jackson 
as an original source of this equation is misleading. The 
monograph [|| does contain derivation of the drift veloc- 
ity. The second term in (l,||4l) in fact is the drift velocity 
Vj^ (Eq.(12.117) from Ref. [g). In the derivation of Vj^ 
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a time-averaging procedure was applied, thus all oscillat- 
ing terms have been canceled. If we are interested only in 
the systematic shift (drift), these terms, of course, do not 
affect the result. However if one wants to calculate the 
radiation of particles, then the oscillating terms cannot 
be ignored. Correspondingly, the velocity v{t) cannot be 
obtained by adding to Vjj a constant vector (of unknown 
direction!) vq, as it is done in Ref. 

It should be noted that in some cases the replacement 
of v{t) by a constant vector should not necessarily lead 
to a wrong results. For example, from the solution of 
the equation for the motion in the constant homogeneous 
field follows that = const, therefore for definition of 
a^, we can ignore the velocity's time-dependence. We 
would like to emphasize that one can make such a state- 
ment only when the solution of the equation motion is 
known. In more complex cases, like the one analyzed 
above, the solution of the equation of motion is a neces- 
sary step. Any a priori assumption that a^(t) = const 
may lead to incorrect results. In this regard, all equa- 
tions in Ref.[Q, which contain are not correct. The 
results of this work can be considered, in the best case, 
as interpolations of two well known cases: <C 
(curvature radiation) and /3j_ (synchrotron radia- 

tion). In these limiting cases the results of Eqs. ( ^ ) and 
( pol) coincide. The reason of this agreement is explained 
above. 

Finally, we note that Ref. |^ contains another incor- 
rectness related to the derivation of Eq.(21,Q). Al- 
though the authors do cite to the classical paper of 
J. Schwinger apparently they have missed to mention 
that not only the final result (i.e. the equation (21, [Q) it- 
self), but also all intermediate steps for derivation of this 
equation (using the same method and the same denota- 
tions) have been already obtained in Schwinger's paper. 



III. INTENSITY AND ENERGY SPECTRUM 
OF RADIATION 

In this section we calculate, within the framework 
of classical electrodynamics, the intensity and the en- 
ergy spectrum of radiation of a charged ultrarelativis- 
tic particle with the perpendicular component of velocity 
P± ^ Pd' If the radius of the curvature of the trajectory 
Rc ^ "fQi the characteristic energy of the emitted photon 
is ~ hc^^ /rQ. The requirement of smallness of this energy 
compared to the energy of the particle mc^7 (the condi- 
tion of applicability of classical electrodynamics) implies 
fi7^ <C mcrQ. We assume that the perpendicular (rela- 
tive to the smooth trajectory) momentum p^ <C Py, but, 
at the same time, p^ ^ mc. Using Eq. (^ , the latter 
inequality can be presented in the form C7^^> uJcVq. This 
gives the following constraint on the Lorentz-factor: 



^c^o /c <C 7^ <C mcro / h . 



(45) 



The condition of smallness of the lower limit compared 
to the upper limit, implies <C mc^, i.e. the mag- 
netic field should not exceed the critical one, Bcr = 
m'^c^/{eh) = 4.4 x 10^^ G. 

The energy radiated away by the particle for a unite 
interval of time, is 



dt~ 3c3 



{V X af/c') = ^ ^ • (46) 



3 m 



Substituting Rc from Eq.(Ea) , we obtain 



(l — 277 COST + 



(47) 



where we introduce a new parameter rj = /3j_//3^. The 
averaged over the period intensity is 

2e^c 7^ 



(48) 



This formula can be written in the form (/) = /curv(l + 
7^^), where /curv is defined in Eq.(|TQ|). It is seen that dur- 
ing the motion in the curved magnetic field the curvature 
radiation has the minimum possible intensity. 

In the considered scenario, we deal with three char- 
acteristic times. These are the time of energy losses 
(cooling time), tcooi ^ mcrQ/(e^7^), the period of "os- 
cillations", T ~ 7/co'c, and the characteristic time of for- 
mation of radiation. At ~ ro/(c7). It is easy to show, 
using Eq.(45) , that At <^ tcooi- Note that this condition 



is always satisfied in the framework of classical electro- 
dynamics. The ratio At/T ~ Uc^q / (c7^) also is small 
compared to 1, as it follows from Eq.diq). Finally, the 
ratio T/tcooi ^ e^7^/(mcco'cro); using the lower and up- 
per limits of the Lorentz-factor from Eq.(p5|), we obtain 



hujc 



a ■ 



^cool 



< a 



mc 

huJr 



(49) 



where a = /hc^ 1/137. 

The smallness of At simplifies significantly the calcu- 
lations of the radiation spectrum. This allows us to ig- 
nore the changes of the particle energy and the curvature 
radius with time, and perform calculations at the fixed 
"prompt" values of these parameters: E{t) = mc^j{t) 
and Rc{t). Then the spectral fiux density of the magne- 
tobremsstrahlung integrated over the all emission angles, 
is described by the well known expression for the syn- 
chrotron regime of radiation (see, e.g. Refs. [^,0) 



P(c^,t) 



2 7r Rc V ^* 



Here 



and 



307^ 



2Rc ' 

F{x)=x K5/s{x')dx', 

J X 



(50) 



(51) 



(52) 



where i^s/s is modified Bessel function. The particle en- 
ergy loss rate, dE/dt = — P{uj^t) duj] the calculation 
of this integral results, as expected, in Eq.(^) 

Since Rc and 7 vary with time, the radiation spectrum 
should be time-dependent as well. Let assume that the 
condition T <C tcooi is satisfied (note that this condition 
does not contradict to Eq.(^). Then, since change of 7 
during the period is small, we can introduce new param- 
eters averaged over time. This in fact has been already 
done at the changeover from Eq.(^ to Eq.(^), assum- 
ing that during the period 7 = const. 

By writing Rc^ in accordance with Eq.(p3[), as 



G{x) 



Rc = ro/q{r], r) , q{r], r) = y^l - 27?cosr + 7^2 ^ (53) 

the averaged spectral flux density can be presented in the 
form 



UJ 



2 IT To V CJ* 



Here 



g(^)='- [\i,,r)F(—^)dr, 
V^*/ ^ Jo \w*(2(r?,T)/ 



(54) 



(55) 



where uj^ = 3c7^/(2ro). G is a function of two variables: 
X = uj/uj^ and r] = ^j_/f3jj. If in Eq. (|55|) we represent the 
function F{x) in the form of Eq.(52), and change the or- 
der of integration, the integral over dr can be calculated 
analytically. This allows us to express the function G{x) 
in the form of a single integral: 



G{x) = X K^/s{x') "^{x/x') dx' , 



where 



1 ix/x') 

- arccos — — - 



1, 



217 



>|i 
< 11 



■ri\ 



(56) 



(57) 



The function G for different values of 77 is shown in 
Fig. ||. At 77 = (the dot-dashed line) we have a nominal 
curvature radiation. The curves are based on numerical 
integration of Eq.(^, using for the function F{x) the 
following analytical approximation from Ref. H, 



F{x) 2.15x^/^(1 + 3.06 x) 

1 + 0.884 + 0.471 x^/^ 
1 + 1.64x2/3 + 0.974x^/3 



xl/6 



(58) 



which provides an accuracy better than 0.2% over the 
entire range of the variable x. In two limits, x ^ 1 and 
X <C 1, the function F{x) can be expressed as 



F(x) ^ 22/^r 



F(x) 



.1/3 



X <Cl , 



X > 1. 



(59) 
(60) 




FIG. 3. The spectral flux density of the small pitch- angle ra- 
diation averaged over time. The thin solid lines (from top to 
bottom) correspond to numerical calculations performed for 
the values of the ratio 77 = I^d — 0? 0-5, 1, 1.5 and 2. The 
dotted lines are calculated using the asymptotic presentation 
of the function G{X) given by Eq.(|6^. The dot-dashed line 
corresponds to the curvature radiation {rj = 0) when the par- 
ticle moves along the soft trajectory. 



Therefore at x ^ 1, the function G(x) differs from F{x) 
only by a factor which does not depend on x: 



1 

x) = - / (l-27?cosr + 7?^)^/^dr xF(x). 



(61) 



In the opposite limit, a; » 1, the integral can be calcu- 
lated by using the standard saddle point method. This 
gives the following asymptotics: 



G: 



2v^ 



■ exp 



1 + 7? 



(62) 



Formally, in this case the saddle point method can be 
applied if x (1 + 7^)^/7?. However, the numerical calcu- 
lations show that in the interval 0.2 < 7/ < 3 the simple 
analytical function given by Eq.(|6^) provides a good ac- 
curacy (better than 10 %) already at x 2 (see Fig. ||). 

Let's consider now the case when the injected accel- 
erated electrons are not strictly unidirectional but are 
distributed within a narrow cone with an opening angle 
<C 1. We assume that all electrons have the same 
energy. For the standard synchrotron radiation, a slight 
change of the pitch angle does not affect the spectrum of 
radiation. However, in the case of the small pitch-angle 
magnetobremsstrahlung, even a tiny spread of pitch- 
angles may result in a significant change of the spectrum. 
It is convenient to express the angular distribution of 
particles through the parameter 7^, which is linked to the 
pitch angle with a simple relation: = r]/3^//3^^. Note 
that at 7^ < 1, the pitch angle 6> <C 1. For calculations 
we should specify the angular distribution of particles. 
Below we consider a Gaussain type distribution: 



g{r]) d7? - e-^'/(2^) dr] , 



(63) 



thus the mean value of the square of the perpendicular 
velocity 

(/?i)=C/?!,. (64) 

For determination of the energy spectrum we have to av- 
erage Eq.([55|) over time r and parameter 77, i.e. calculate 
the following integral 

(OM)=p,»(,)£f,(,,r,F(^), (65) 

The energy spectra calculated for several values of ( 
are shown in Fig. ^. Apparently, the case of ^ = corre- 
sponds to the curvature radiation. 

As above, using the saddle point method, we can find 
the asymptotics for (G) at large x. Rather cumbersome 
calculations, which we omit here, leads to the following 
result 

where 

u^^ C^/^x^'^ - C-2/3a;i/3 + i. . (67) 

This equation is derived under the following conditions: 
X ^ and x ^ !/(". However, the comparison of Eq.(pj) 
with accurate numerical calculations show (see Fig.B), 
that for C ~ 1 this analytical presentation gives correct 
results already at x > 0.5. 

Although the assumed Gaussian type angular distri- 
bution of electrons seems to be a quite reasonable and 
natural choice in the considered scenario, it would be 
interesting to investigate the dependence of the radia- 
tion spectrum on the specific angular distribution of the 
electron beam. In particular in Fig.|^ we show the radi- 
ation spectra calculated for the uniform distribution of 
electrons within fixed opening angles of the beam. The 
results shown in Fig.|| (solid curves) are obtained for the 
same values of the mean square of the perpendicular ve- 
locity as in Fig. ^. In Fig.|| we show also the asymptotic 
solutions given by Eq.(p6|) for the spectra calculated for 
the Gaussian type angular distribution (the same dotted 
curves shown in Fig. ^ . The comparison of the results 
shown in Figs. ^ and || indicates that the value of the mean 
square of the perpendicular velocity generally character- 
izes the energy spectrum of radiation independent of the 
details of the specific small pitch-angle distribution of 
electrons. 

Note that the spectra of the small-pitch angle radiation 
do not contain the characteristic exponential term e~^, 
as in the case of the synchrotron or curvature radiation, 
but show a harder behavior. This can be seen from the 
asymptotics Eq. (|67|) . To demonstrate the tendency of 
steepening of the energy spectrum beyond the maximum 
at X ~ 1, in the table below we show the evolution of 
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{G{x)) 

z 1 1 — I — I I I I I I 1 1 — I — I I I I I I 1 1 — I — I I I I I I 




FIG. 4. The spectral flux density of the small-pitch angle 
radiation averaged over the period and the pitch-angles of 
particles, assuming a Gaussian type of distribution given by 
Eq.di^). The solid lines correspond to different values of ^ = 
{P±) / Pd' I5 3 and 5 (from bottom to top). The dotted lines 
are calculated using the asymptotic analytical presentation 
given by Eq.(^). The dash-dotted line corresponds to the 
curvature radiation (^ = 0). 



the local slope of the spectral fiux density ( the so-called 
spectral index a) which is determined as 

a = -d{\og{G))/d{\oguj) . (68) 

It is seen that while in the case of curvature radiation 
(C = 0), the spectrum becomes very steep already at 
X > 3, in the small-pitch angle radiation regime with 
(^ = 2 or (^ = 3, a relatively hard spectrum can extend 
up to X ~ 10. 



>^ 


2.0 


3.0 


5.0 


10. 


20. 


30. 





1.65 


2.62 


4.59 


9.56 


19.5 


29.5 


0.5 


1.17 


1.67 


2.56 


4.45 


7.56 


10.2 


1.0 


0.96 


1.37 


2.10 


3.63 


6.14 


8.28 


2.0 


0.75 


1.08 


1.67 


2.91 


4.93 


6.65 


3.0 


0.64 


0.93 


1.45 


2.54 


4.31 


5.82 



IV. DISCUSSION 



In the context of the motion of an ultrarelativistic 
charged particle in the strong magnetic field with curved 
field lines, in the vicinity of the latter always exist tra- 
jectories without fast oscillations of the curvature radius. 
We call them smooth trajectories. The angle between the 
smooth trajectory and the field line appears equal to Pjj 
and is determined by Eq.(^. For the conditions dis- 
cussed in this paper, Pjj <C 1. This trajectory is distinct 
in the sense that, strictly speaking, the curvature radia- 
tion takes place when the charged particle moves along 
this trajectory, but not along the magnetic field line as 
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<C(x)> 



1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 




^^^^^^^ ' ' ^^^^^^<^ 






















1 





FIG. 5. The spectral flux density of the small-pitch angle ra- 
diation averaged over the period and the pitch-angles of par- 
ticles, assuming a uniform angular distribution of electrons 
within the opening angle of the beam. The solid lines corre- 
spond to different values of the mean square of the perpen- 
dicular velocity C = 1? 3 and 5 (from bottom to 
top). The dotted lines are the same as in Fi^. They are 
calculated for the Gaussian type of angular distribution us- 
ing the asymptotic analytical presentation given by Eq.(|66|). 
The dash-dotted line corresponds to the curvature radiation 

(C = 0). 



it is broadly accepted and discussed in the literature. 
In fact, if the velocity vector is initially parallel to the 
magnetic field line, no radiation is possible hy definitionl 
The curvature radiation appears when the velocity vector 
becomes parallel to the smooth trajectory. On the other 
hand, since the curvature radius of the smooth trajectory 
Rc is very close to the radius of the field curvature ro, 
the intensity and spectral distribution of the curvature 
radiation are nevertheless correctly described by the for- 
malism applied in the original papers (see e.g. B) where 
the curvature radiation has been linked to the motion of 
particles along the magnetic field lines. 

Another common misconception in the literature is 
connected with the belief that the radiation in the strong 
magnetic field results in prompt damping of the perpen- 
dicular component of motion. However, if the perpendic- 
ular momentum is ultrarelativistic, the parallel and per- 
pendicular components of the momentum decrease (due 
to the radiation) by the same law. A strict proof of this 
statement is possible in the case of homogeneous mag- 
netic field. However, in the context of consideration of 
direction of the recoil momentum, this feature can be 
extended to the curved magnetic field for the particle 
motion along the trajectories close to the smooth tra- 
jectory. The fast damping of the perpendicular motion 
occurs only if the perpendicular component of momen- 
tum is non-relativistic. This circumstance might have 
important practical implications for radiation of highly 
magnetized relativistic objects. 

If initially the velocity vector of the particle is not di- 



rected strictly along the smooth trajectory, but at a small 
angle {0 <C 1), then for ^7 1 the perpendicular mo- 
mentum is ultrarelativistic and the particle moves along 
the line spiraHng around the smooth trajectory (but not 
the field line). In particular, this happens when ini- 
tially the velocity is parallel to the field line, provided 
that jPjj ^ 1. Then, at the initial moment of time 
a = and Rc = 00. The important difference from the 
motion in the homogeneous magnetic field is the time- 
dependence of the curvature radius and the square of 
acceleration. The radiation in this regime is principally 
different from both the synchrotron and the curvature 
radiation regimes. This new regime of the magneto- 
bremsstrahlung, which surprisingly has been missed in 
the previous studies, would be appropriate to call small 
pitch-angle radiation. 

The averaged over the period intensity of radiation is 
given by Eq .(p8|) . It can be presented in the standard 
form of Eq.([46|), if we introduce an effective curvature 
radius i?eff = ^o/vI + t^. However, if we want to have 
a standard exponential term in an explicit form (as 
in the case of the synchrotron or curvature radiation) in 
the asy mptotic presentation for the spectrum given by 
Eq.(p2|), the effective curvature radius should be defined 
differently, namely i?eff = ro/{l -\- r]). This is another 
indication of the difference between the small-pitch angle 
radiation from the synchrotron and curvature regimes of 
radiation. 

The characteristic feature of the small-pitch angle ra- 
diation is its strong dependence on the particle angular 
distribution. Even in the case of a very narrow angular 
distribution, the radiation spectrum could differ signif- 
icantly from the unidirectional beam. This is demon- 
strated in Fig.^ and || where the radiation spectra aver- 
aged over the pitch-angles of particles, assuming Gaus- 
sian type and uniform angular distributions of particles 
are presented. For the Gaussian type angular distribu- 
tion, the spectral flux density can be presented in a sim- 
ple analytical form by Eq.(|6q) which provides a good 
accuracy already at x > 0.5. It is seen that the small 
pitch-angle radiation predicts significantly harder energy 
spectra compared to the curvature radiation ((" = 0). 

This interesting feature might be a key for interpreta- 
tion of the recent observations of pulsars by Fermi LAT 
which indicate that the energy spectra of some pulsars, 
in particular the Crab pulsar [|[ , in the cutoff region are 
significantly harder than as predicted by the cur- 
vature radiation models. The energy spectrum of the 
pulsed emission of the Crab reported in Ref. p can be 
readily described by Eq.(^) assuming C = 1 and a quite 
reasonable value of the ratio — 10^^ cm~^. 

We should note, however, that the formalism devel- 
oped in this paper cannot be directly used for inter- 
pretation of experimental data. It can be applied only 
at the initial stages, when the energy losses of particles 
are negligible. Otherwise, for a correct treatment of the 
time-evolution of particles and their time-dependent ra- 
diation, one should solve kinetic equations that describe 
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FIG. 6. The particle trajectories on the planes (x, z) (left 
panel) and {x,y) (right panel) for angles 19 = 10° and 15°. 
The motion is infinite, the direction of the final velocity de- 
pends on the initial angle 19. The circle at the origin of coor- 
dinates is the surface of the "star" . 




the energy and space distributions of particles taking into 
account the particle energy losses, changes of their tra- 
jectories, etc. This issue we will discuss elsewhere. 

Finally, we note that the results of this work are ob- 
tained within the framework of classical electrodynamics 
when the quantum effects can be ignored. If the condi- 
tion of Eq.(45) is violated, then calculations should be 
conducted using the methods of quantum theory. If in 
the initial state the electrons occupies the first level of 
Landau, then the standard method of calculations used 
in the quantum theory of synchrotron radiation (see e.g. 
Ref. ||l^) is not applicable since the perpendicular mo- 
tion cannot be treated as quasi-classical. 



Appendix A: Trajectories in the field of magnetic 
dipole 

Here we present for illustration the results of numeri- 
cal calculations of trajectories in the field of the magnetic 
dipole. The initial Cartesian coordinates of the particle 
(in units of R^) we set on the sphere of radius in 
the form r(0) = (sini^, 0, cosi9); the initial velocity is 
directed along the field hue, which starts from the point 
r(0). It is convenient to introduce, instead of t, a dimen- 
sionless parameter r = t eBo/{ncj)^ where Bq is the field 
in the initial point of the trajectory, and write down the 
equation of motion in the following form 



FIG. 7. The same as in Fig. fel but for angle = 20°. 



dr 

d^ 



(Al) 



Here u = (5/13 is the unite vector in the direction of 
velocity, b = B{r)/ Bq^ 



cBqR^ 



(A2) 




FIG. 8. The same as in Fig. |, but for angles i9 = 30° 
45°. 



and 



Note that only one parameter, ^, enters into this equa- 
tion. In accordance with the condition of Eq.(|8|), the 
values of 1 — P and ^ are small compared to unity. We 
adopted ^ = 10~^, and performed the calculations with- 
out taking into account the effect of radiative friction. 

If initially we put the particle on the "pole" and direct 
the velocity towards the axis the particle apparently 
will move uniformly along that axis. If at t = the 
particle is located close to the "pole" , the motion will be 
not rectilinear, but it will remain infinite (see Fig. ||). 

With an increase of i9 the character of motion is dra- 
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matically changed. For example, at = 20° the particle 
is trapped by the magnetic field (see Fig. 0). There is 
no sharp transition between these two cases. At some 
angles close to i9 = 20°, the particle may make several 
turns around the star, with a quite chaotic features of 
motion, and then go to infinity. 

At large angles the trajectories are again changed 
quantitatively, and appear close to the field lines (see 
Fig. H). In all cases the initial segment of trajectory ap- 



pears close to the field lines; the acceleration on that 
segment is shown in Fig. |[ Therefore, for the chosen 
initial conditions the curvature radius of the trajectory 
significantly differs from the curvature radius of the field 
line. The initial segments of the trajectory are helical- 
like lines, but at ^ ~ 10~^ the shift in the perpendicular 
direction appears of the order of 10~^. Therefore within 
this uncertainty, the trajectories of particles shown in 
Figs. 1^ - H coincide with the smooth trajectories. 
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